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We consider the J/\t suppression at a high energy heavy ion collision. An ideal gas of massive 
hadrons in thermal and chemical equilibrium is formed in the central region. The finite-size 
gas expands longitudinally in accordance with Bjorken law. The transverse expansion in a 
form of the rarefaction wave is taken into account. We show that J/ty suppression in such an 
environment, when combined with the disintegration in nuclear matter, gives correct evaluation 
of NA38 and NA50 data in a broad range of initial energy densities. 
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I. INTRODUCTION 

Since the paper of Matsui and Satz pi] there is a steady interest in the problem of J/^ suppression in a 
heavy-ion collision. The question is if this suppression can be treated as a signature for a quark-gluon plasma 
(QGP) or if it can be explained by J/vE" absorption in a hadron gas which appears instead of the QGP in the 
central rapidity region (CRR) of the collision In the following paper we shall explore the latter possibility. 

The existence of the QGP — the novel phase of hadron matter, has been predicted upon lattice QCD calcu- 
lations (for a review see |Q and references therein). One gets from those calculations the critical temperature T c 
in the range of 150 - 270 MeV. The upper limit belongs to a pure SU (3) theory, whereas adding quarks causes 
lowering of T c even to about 150 MeV (Nf > 3). All lattice calculations have been performed for zero quark and 
baryon chemical potential only. Additionally to the fact that lattice estimations could not be treated literally, 
in relativistic heavy ion collisions there is a finite baryon chemical potential. As a result, the above-mentioned 
evaluations of T c can be understood only as a measure of the order of the magnitude of real T c . Therefore, the 
assumption that at initial temperatures around 200 MeV (as in our model) a hadron gas still exists is realistic 
entirely. As far as the hadron gas is considered itself in that range of temperature, it was shown in || that the 
temperature increases with energy density when continuum excitations (string degrees of freedom) are not taken 
into account. Continuum excitations gives the limiting temperature. This reproduces results of the Hagedorn 
bootstrap model H . The thermodynamical analysis can also be based on particle ratios [110] and this gives the 



temperature as in our model. The same result is obtained from the microscopic cascade model (see e.g. 11 ). 

In the following paper we shall continue our previous investigatio ns |12[ ] into the problem of J/*ff suppression 
observed in a heavy-ion collision (for experimental data see e.g. |l3| , |l4 | and references therein). Now, we shall 
focus on the dependence of the suppression on the initial energy density reached in the CRR. 

In our model, J/^f suppression is the result of a cc state absorption in a dense hadronic matter through 
interactions of the type 

cc+h — >D + D + X, (1) 

where h denotes a hadron, I? is a charm meson and X means a particle which is necessary to conserve the 
charge, baryon number or strangeness. 

First, the J/\I> is absorbed by nucleons in colliding nuclei [ ^5[ . Later, the J/ty is absorbed by secondary 
hadrons produced in the collision. In the simpler scenario those secondary hadrons expand longitudinally along 
collision axis Jl^Jl7]] and the time evolution is given by Bjorken's scaling dynamics fl8f| . The importance of 
secondaries scattering becomes more and more important at higher energies since secondary production grows 
with energy as well. We assume that at very high energies secondaries form a dense hadronic gas in a state 
of thermodynamical equilibrium. A chemistry will be given by equation of state of an ideal gas consisting of 
different species of massive hadrons. 

A problem of thermal and chemical equilibration of an hadronic fireball formatted in a heavy collision is still 
far from being solved. Recent results lTl|J^,|l9|, based on parton cascade models, suggest that in big systems the 
equilibration time for heavier particles is longer than for lighter particles. The thermal equilibrium is established 
quickly, in about 5 fm/c, much faster than the chemical equilibrium. 
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To keep our model as simple as possible we consider the one temperature model. This allows us to reduce the 
number of external parameters. It is obvious that any new parameter would result in a better fit to experimental 
data. 

The hadronic matter is in a state of an ideal gas of massive hadrons in thermal and chemical equilibrium 
and consists of all species up to fl~ baryon. Time evolution is given here by conservation laws combined with 
assumptions about the space-time structure of the system. A corresponding equation of state of the ideal gas 
makes then possible to express gas parameters such as temperature and chemical potentials as functions of time. 

An ideal gas of real hadrons has a very interesting feature: it cools much slower than a pion gas when 
expands longitudinally. We have checked numerically that for the initial energy densities eo corresponding to 
initial temperatures To of the order of 200 MeV and for the freeze-out T/. Q . not lower than about 100 MeV, 
the time dependence of the temperature of the expanding gas still keeps the well-known form T(t) = T ■ t~ a 
(we put to = 1 fm). Only the exponent a changes from ^ for massless pions to the values 5^-5^3 for massive 
realistic hadrons. As a result, the time of the freeze-out tf. is much greater for the hadron gas than for the 
pion one. For instance, when we take To = 200 MeV and T/. . = 140 MeV we obtain tf. a . — 7.37 fm (a = J^) 
for the hadron gas and tf. a . — 2.9 fm (a = -|) for the pion gas. The lower I/. ., the stronger difference. For 
T = 200 MeV and T f . . = 100 MeV we have t f . . = 48.5 fm (a = -^) for the hadron gas and t f . . = 8 fm 
(a = i) for the pion gas. This has a direct consequence for J/^S suppression: the longer the system lasts, the 
deeper suppression causes (see Fig. |l|) . 
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FIG. 1. Comparison of suppression of the pure J/\t's in the cooling hadron gas for two values of the power a in 
approximation T(t) = To ■ t~ a : a = \ (solid) and a = (dashed). 

We are going to calculate a survival factor for J '/vff when new, more realistic conditions are taken into account. 
We consider a hadronic gas which is produced in the CRR region. This gas expands both longitudinally and 
transversely. The longitudinal expansion is a traditional adiabatic hydrodynamical evolution p8[ , the transverse 
expansion is considered as the rarefaction wave. An initial energy density depends now on the impact parameter 
b and on the geometry of the collision. 

The concept of the description of J/\l/ absorption in the hot hadronic matter is based on the application 
of the relativistic kinetic equation as it was postulated in ||. As the hadronic medium a pion gas was taken 
into account there. Our hadronic medium consists of hadrons from the lowest up to f2 _ as constituents of the 
matter in the CRR. This gives different physical properties important for the absorption processes as compared 
to the absorption in a pion gas — cooling time, threshold effects, etc. It allows us to consider here also a case 
with non-zero baryon number density. We take into account here some physical effects which were neglected 
or not fully treated in (^]. Thresholds for J/ty — hadron reactions as well as relative velocities are examined 
completely, but both these effects are ignored in final calculations in In addition to the CRR also J/^ 
absorption in the nuclear matter is considered simultaneously here. In the final estimations of J survival 
fraction, the Woods-Saxon nuclear matter density distribution is used in the presented model. 

As far as the hydrodynamic is concerned, the full solution of hydrodynamic equations (for the baryon number 
density equal to zero) with the use of the method developed in p0[ was explored in J5j , but with different initial 
temperature profile. Here, for simplicity of numerical calculations we assumed the uniform initial temperature 
distribution with the sharp edge at the border established by nuclei radii (see Fig.0). For such an initial 
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distribution, and for a central collision and the baryon number density equal to zero, hydrodynamic equations 
are solved numerically in [ pp| . The resulting evolution can be decomposed into the longitudinal expansion inside 
a slice bordered by the front of the rarefaction wave and the transverse expansion which is superimposed outside 
of the wave. But the temperature decreases very quickly outside the wave (at least at times not greater than the 
half of the freeze-out time), so it is reasonable to ignore the influence of the flow outside of the rarefaction wave 
as we do in our model. Since we are dealing with small but nevertheless non-zero baryon number densities, we 
assume that the full hydrodynamic evolution looks like the same in this case. For the longitudinal expansion — 
the most important part of the hydrodynamical evolution, this assumption is confirmed by our investigations 
into the time dependence of the temperature of the hadron gas, where the same pattern (see (^2|)) has been 
obtained approximately as in the case of a massless baryonless ideal gas |2l],|22j . 

The dependence of J/^f suppression on the impact parameter b is treated here more rigorously than in |||. 
In fact the only approximation of the survival factor (with respect to its exact expression for the case of the 
hnite-size effects and the transverse expansion taken into account) is done by treating the time J/^P needs to 
escape from the hadron matter in the transverse plane on the average (over the initial J/ty positions). The 
effective freeze-out time tf (denoted by t final herein) is therefore the function of both the impact parameter and 
J/ty transverse momentum, whereas tf is evaluated only for given b (for some assumed average J/^ transverse 
momentum, namely equal to 1 GeV) in 



II. THE EXPANDING HADRON GAS 

For an ideal hadron gas in thermal and chemical equilibrium, which consists of I species of particles, energy 
density e, baryon number density ub, strangeness density n$ and entropy density s read (h — c — 1 always) 



2lr2 i=i J " <<M> 



£ =^>.(2 Sl + l) / dp ^— , (2a) 



n B 



n s 
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£(2 Si + l)/ dp , (2b) 

2n ti J° exp{^}+ 5l 

£(2^ + 1)/ dp , (2c) 

2ir "tt J o exp{^}+ ffj 
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where Ei = {mf + p 2 ) 1 / 2 and mi, Bi, Si, fa, Si and gi are the mass, baryon number, strangeness, chemical 
potential, spin and a statistical factor of specie i respectively (we treat an antiparticle as a different specie). 
And fa = Bi/is + Sifis, where jib and are overall baryon number and strangeness chemical potentials 
respectively. 

We shall work here within the usual timetable of the events in the CRR of a given ion collision (for more 
details see e.g. ||). We fix t = at the moment of the maximal overlap of the nuclei. After half of the time the 
nuclei need to cross each other, matter appears in the CRR. We assume that soon thereafter matter thermalizes 
and this moment, to, is estimated at about 1 fm |i 18|. Then matter starts to expand and cool and after 



reaching the freeze-out temperature it is no longer a thermodynamical system. Wc denote this moment as tf. . 
As we have already mentioned in the introduction, this matter is the hadron gas, which consists of all hadrons 
up to Q~ baryon. The expansion proceeds according to the relativistic hydrodynamic equations and for the 
longitudinal component we have the following solution (for details see e.g. ]l8| , p3| ]) 

S W = — , n B{t) = , (3) 
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where so and n° B are initial densities of the entropy and the baryon number respectively. The superimposed 
transverse expansion has the form of the rarefaction wave moving radially inward with a sound velocity c s and 
the transverse flow which starts outside the wave [jl8],^0| . But as it has been mentioned in the introduction, since 
the temperature decreases rapidly outside the wave (for the most important times for absorption processes), we 
shall ignore possible J/\& scattering there. It means that J/vf suppression is slightly underestimated here. 

To obtain the time dependence of temperature and baryon number and strangeness chemical potentials one 
has to solve numerically equations ( pb| - pc| ) with s, ns and n$ given as time dependent quantities. For s(t), 
n-B(t) we have expressions ([|) and n$ = since we put the overall strangeness equal to zero during all the 
evolution (for more details see 

The sound velocity squared is given by c 2 a = ^ and can be evaluated numerically |2l],^2| . 



III. J/* ABSORPTION IN HADRONIC MATTER 



In a high energy heavy-ion collision, charmonium states are produced mainly through gluon fusion and it 
takes place during the overlap of colliding nuclei. For the purpose of our model, we shall assume that all cc pairs 
are created at the moment t — 0. Before the fusion, gluons can suffer multiple elastic scattering on nucleons 
and gain some additional transverse momentum in this way [p4|-p6[ . This manifests for instance in the observed 
broadening of the pt distribution of J/4" |27|]. Following ^8|, we express this effect by the transverse momentum 
distribution of the charmonium states 



5(PT,eo) 



2 P7 



(Pl)^(eo) 



exp 



Pt 



(Pl)^(eo) 



(4) 



where (pj) j^Xeo) is the mean squared transverse momentum of J/* gained in an A-B collision with the initial 
energy density eg. The momentum can be expressed as (for details see ]29|]) 



(Pt)jM^ = ( Vt) P j% + 



K-e 



(5) 



with K = 0.27 fm 3 GeV and (pt) P j% 



1.24 GeV 2 taken from a fit to the J/* data of NA38 Collaboration 
is normalized to unity and is treated as the initial momentum distribution of charmonium 



The expression in 
states here. 

For the simplicity of our model, we shall assume that all charmonium states are completely formed and can be 
absorbed by the constituents of a surrounding medium from the moment of creation. It means that we neglect 
a whole complex process of J/\& formation as presented in |30| , |3l|| . The main feature of the above-mentioned 
process is that, soon after the moment of production, the cc pair binds a soft gluon and creates a pre-resonance 
cc — g state, from which, after a time of the order of 0.3 fm, a physical charmonium state is formed. This means 
that the possible nuclear absorption of charmonium is, in fact, the absorption of the cc — g state. But the 
latest has the cross-section a a i, s = 7.3 mb, which is much higher than J/ty — Nucleon absorption cross-section 
&il>N =3 — 5 mb obtained from p-A data |]l5| , p7| , j32"[ . This justifies our assumption: taking into account cc — g 
absorption instead of charmonium disintegration in the nuclear matter would only strengthen J/\E' suppression. 

According to the above assumption, charmonium states can be absorbed first in the nuclear matter and soon 
later, when the matter appears in the CRR, in the hadron gas. Since these two processes are separate in time, 
J/vP survival factor for a heavy-ion collision with the initial energy density eo, m &y be written in the form 



Af(e ) = jV„. m .(e ) • Afh. g .(e ) 



(0) 



where JV n , m .(^o) and Nh.g.( e o) are J/^ survival factors in the nuclear matter and the hadron gas, respectively. 
For M n .m.( e o) w e have the usual approximation J3^—|35[] 

Nn.m.(eo) — exp{-cr^Ar/9oi} , (7) 

where po is the nuclear matter density and L the mean path length of the J/^F through the colliding nuclei. For 
the last quantity, we use the expression given in p5[: 



PoL(b) 
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■T A (s) + ——T B (s-b) 
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(8) 
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where T^g(b) = J d 2 s Ta(s)Tb{s — b), T^(s) = j dzp^s, z) is the nuclear density profile function, p^(s,z) 
the nuclear matter density distribution (normalized to A) and b the impact parameter. How to obtain eo as a 
function of b will be presented further. 

To estimate Afh.g.( e o) we follow the idea presented in ||, but now generalized to the case of the gas which 
consists of different species of particles. We shall focus on the plane z = (z is a collision axis) and put J/^f 
longitudinal momentum equal to zero. Now the p^-dependent J/^ survival factor Nh.gXPT) is given by (for 
details see [Q) 

Nh.gXPT) = J ^ 2 s7o(s,p T )exp j- dtY^J -^^h{q,t)^Vrel,J^^ , (9) 

where the sum in the power is over all taken species of scatters (hadrons), p" — (E,pr) and q\ = {E[,q) are 
four momenta of J/$ and hadron specie i respectively, v = pr/E is the velocity of the former, Ui states for the 
absorption cross-section of J/^S — hi scattering and v r ei,i is the relative velocity of hi hadron with respect to 
J/\P. When M and rrii denote J/ty and hi masses, respectively [M — 3097 MeV), u re z,t reads 

/ m 2 M 2 V 

The upper limit of the time integral in (Q) , tf, is equal to tf. , or to t esc - the moment of leaving by a given 
J /'fy of the hadron medium, if the final-size effects are considered and t esc < tf . For cr^ we have assumed that 
it equals zero for (jf -\-q" ) 2 < (2tod + mx) 2 and is constant elsewhere (to_d is a charm meson mass, mjj = 1867 
MeV). For hadron specie i we have usual Bose- Einstein or Fermi-Dirac distribution (we neglect any possible 
spatial dependence here) 

mt) = f&,t)= + ; — . (id 

eXP \ T(t) J+gi 

In the following, we shall consider only J/^S initial distribution /o(s,pr) that factorizes into fo(s)g(pr) and 
the momentum distribution g(pr) will be given by (^|) . We assume at the first step that the transverse size of 
the hadron medium is much greater than tf. Q , and also much greater than the size of the area where fo(s) is 
non-zero. Additionally we assume that fo(s) is uniform and normalized to unity. Note that the first assumption 
overestimates the suppression but the second, in the presence of the first, has no any calculable effect here. As 
a result, Mh.g.ipr) simplifies to 

{f^f-°- ^ r d^Q p 1 

To obtain Afh.g.( e o) ° n e needs only to integrate (|l2|) over px- 



Nh.g.{ e o)= I dp T g(p T ,e ) -exp I - I dt V / JrLfi($, t)(XiV re i,i^% \ ■ (13) 



(2tt) 3J,w ' EEi 



Now we would like to take the final-size effects and the transverse expansion into account in our model. To 
do this directly, we would have to come back to the formula given by (||) and integrate it, instead of (|lj) , over 
Pt- But this would involve a five-dimensional integral (the three-dimensional integral over q simplifies to the 
one-dimensional one, in fact) instead of the three-dimensional integral of (|l3|). Therefore, we need to simplify 
in some way the direct method just mentioned above. We shall define an average time of leaving the hadron 
medium by J/'S's with the velocity v produced in an A-B collision at impact parameter 6, (t esc )(b,v). Then, if 
this quantity is less than tf. a ,, we will put it instead of tf. a . as the upper limit of the integral over t in (|l3|). 
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FIG. 2. View of a Pb-Pb collision at impact parameter b in the transverse plane (z = 0). The region where the nuclei 
overlap has been hatched and its area equals S e ff. 

Let us consider an A-B collision at impact parameter b. Since we will compare final results with the latest 
data of NA50 which are for Pb-Pb collisions @, we focus on the case of A=B here. So, for the collision at 
impact parameter b we have the situation in the plane z = as presented in Fig.^, where S e ff means the area 
of the overlap of the colliding nuclei. We shall assume here, that the hadron medium, which appears in the 
space between the nuclei after they crossed each other also has the shape of S e ff at to in the plane z — 0. And 
additionally, the transverse expansion starts in the form of the rarefaction wave moving inward S e ff at to- Then, 
for a J/^l which is at r £ S e ff at the moment to and has the velocity v we denote by t esc the moment of crossing 
the border of the hadron gas. It means that t esc is a solution of the equation | d + v(t — t ) |= R A — c s (t — to), 
where R A — ?~o ■ A s is the nucleus radius and d = r — b for the angel between r and v such that the J/\& will 
cross this part of the edge of the area of the hadron gas which was created by the projectile and d = r in the 
opposite. Having obtained t esc , we average it over the angel between r and v, i.e. we integrate t esc over this 
angel and divide by 2ir. Then we average the result over S e ff with the weight given by 

,-f. T A {r)T B {r- b) 
PjMr) = T AB (b) (W) 

and we obtain (t esc )(b,v). So, the final expression for Afh.g.( e o) when the transverse expansion is taken into 
account reads 

A4. 3 .(eo) = J dp T g(pr,eo) • exp j - ^ J ^^^^ aiVrel '^^i\ ' ( 15 ) 

where t fina i = min{(t esc ),tf, ,}. 



IV. THE ENERGY DENSITY IN THE CRR 

We compare our theoretical estimations for J/^ survival factor with the experimental data JLJ] presented as 
a function of eo- Usually, this quantity is estimated from the well-known Bjorkcn formula 

60 = a q — T ' 16 

Arjbeffto 

where A77 is the pseudo-rapidity range and Et is the neutral transverse energy. 

In further considerations we will need the formula for the number of participating nucleons as a function of 
impact parameter b, which is given by the rough approximation (commonly used in the early nineties) 

N part (b)= [ d 2 s\T A (s) + T B {s-b)\ , (17) 
or in the term of the number of "wounded" nucleons J30| , p5t 



J/ty suppression in colliding nuclei: statistical model analysis 



7 



-V„ ,/(/'!= |^T,(?i[l- l-^-T B (s-b) 



d z sT B {s-b)\l- l-^-T A (s) 



(18) 



The both expressions are depicted in Fig.|^. Note that N part (b) is estimated for the uniform nuclear matter 
density with tq = 1.2 fm and tq = 1.12 fm, whereas N woun d(b) is evaluated for the Woods-Saxon nuclear matter 

With the use of N woun d(b), the relation between Et and 



density distribution with parameters taken from 
b is established @|| 

Et = 



id{b) , 



(19) 



where q = 0.4 GeV. We could put ( |19| ) into (|16|) to obtain the dependence between e and b but the ratio 
N S°Jf(b) ^ * s divergent when 6 — > -fiU + i?s (i?A and i?s are radii of a projectile and a target, respectively). 
From Fig. || we can see that in the case of Pb-Pb collisions, N woun d(b) do not differ substantially from N part (b) 
with ro = 1.2 fm (besides the low b region). Therefore we can assume that for Pb-Pb collisions of NA50 the 
following approximation is valid: 



Having put rt2Q) into ( 



Et — q ■ Np ar t 
we obtain eo as a function of b 



eo(b) 



Npqrt(b) 
Seff(b) 



(20) 



(21) 



where we have also used the value An = 1.2 of NA50 H]. The above function is depicted in Fig. [I], together 
with eo(b) obtained from eq{Et) with the use of (p^). The dependence of eo on Et has been extracted directly 
from NA50 data Q. 
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FIG. 3. Number of participating nucleons as 
a function of b for a Pb-Pb collision, estimated 
as: N P art(b) for the uniform nuclear matter den- 
sity and ro = 1.2 fm (solid), ro = 1.12 fm 
(dashed); N woun d(b) for the Woods-Saxon distribu- 
tion (short-dashed). 
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FIG. 4. The initial energy density eo in the CRR 
for Pb-Pb collisions as the function of impact param- 
eter 6 extracted from the NA50 data Q (solid) and 
obtained from ( pl| ) (short-dashed). 



We shall later see that the low b region (central collisions) is crucial for the understanding of experimental 
data. Also different plots in Fig. [I| will lead to different ranges of initial energy densities obtained from formula 
(§3) or given by NA50 Collaboration gl. 



V. RESULTS 



To evaluate formula e (|l3|) and ( [l5|) we have to know T(t), /is(i) and ns(t) and how to obtain these functions 
was explained in Sect. |II|. But to follow all that procedure we need initial values sq and n B . To estimate initial 
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baryon number density n B we can use experimental results for S-S p7| or Au-Au |l^,[38| collisions. In the 
first approximation we can assume that the baryon multiplicity per unit rapidity in the CRR is proportional 
to the number of participating nucleons. For a sulphur-sulphur collision we have dNs /dy = 6 |57| and 64 
participating nucleons. For the central collision of lead nuclei we can estimate the number of participating 
nucleons at 2 A = 416, so we have dNs/dy = 39. Having taken the initial volume in the CRR equal to irR\ ■ 1 
fm, we arrive at n B = 0.25 fm~ 3 . This is some underestimation because the S-S collisions were at a beam 
energy of 200 GeV/nucleon, but Pb-Pb at 158 GeV/nucleon. From the Au-Au data extrapolation one can 
estimate n° B = 0.65 fm -3 JlQ] . These values are for central collisions, and for the higher impact parameter (a 
more peripheral collision) the initial baryon number density should be much lower. So, to simplify numerical 
calculations we will keep n B constant over the all range of b and additionally, to check the possible dependence 
on n B , we will do our estimations for n B substantially lower, i.e. n° B = 0.05 fm -3 . 

Now, to find sq, first we have to solve (Ea - Ed) with respect to T, ps and /j,b, where we put e — €q, ub = n % 



and ns — 0. Then, having put T, ps and Pb into (2d) we obtain sq. Finally, expressing left sides of (2b,2q) 
by (^) and after then solving 



2d) numerically we can obtain T, /is and /is as functions of time. In fact, 
evaluating formulae ( |l3| ) and (^[) we do the following: first, we calculate T — T(t) which turns out to be very 
well approximated by the expression 

T(t) = T • r° (22) 

and then we put this approximation into (|l3|) and ( p"5| ) . And for ps(t) and PB(t) in ,fi(q,t) we put solutions 
of ( pb| , |2c[ ) where is given by (||), T by (|22|) and n$ — 0. But the exponent a in ( |22| ) has proven not to be 
unique for the whole range of To considered here. One gets different values of the initial energy density eo for 
different values of the impact parameter b and for different geometry of the collision process. So b dependent 
a gives also b dependent freeze-out time The density ep is extracted from the dependence represented by 
the solid line in Fig.|| for different values of b and Eqs. (|2a| - ^cj) are solved. We have evaluated the suppression 
factor up to eo = 3.7 GeV/fm 3 . This gives the maximal possible To, To jnax , equal to 221.8 MeV (for n B = 0.65 
mi" 3 ), 226 MeV (for n° B = 0.25 fm" 3 ) or 226.7 MeV (for n° B = 0.05 far 3 ). 

This procedure allows to evaluate J/vf survival factor given by ( |l3| ) . Because of the lack of data, we shall 
assume only two types of the cross-section, the first, erf,, for J/^-baryon scattering and the second, a m , for 
J/^P-meson scattering. For cr;, we put o\, = ctj/^n- As far as a m is concerned, we assume that this cross-section 
is 2/3 of the corresponding cross section for baryons, which is due to the quark counting. In the following, we 
will use values of J/vE' — Nucleon absorption cross-section aj/^,N = 3 — 5 mb obtained from p-A data |l5|f27| , |3"2| ] . 
At the beginning, to illustrate how the value of power a influences J/^P suppression we present in FigTjlf two 
results: the first for a = | (what is the exact value for a free massless gas) and the second for a = (what is 
the approximate value for the hadron gas and To = 200 MeV). We can see that the suppression improves more 
than twice for the highest eo indeed. 

To make our investigations more realistic we have to take into account that only about 60% of J/^ measured 
are directly produced during collision. The rest is the result of \ (~ 30%) and ip' (~ 10%) decay [BlJ. Therefore 
the realistic J /iff survival factor should read 

Af(eo) = 0.6Mj/4e ) + 0.3A^(e ) + 0.1AVM , (23) 

where Afj/ip(eo), Af x (co) and Af^'(eo) are given also by formulae (^[jl5|) but with (Pr)j^(e) = 
(Pt) j/*( e )' (Pt>x S ( £ )' (PT)$?( e )> K J/-4> = K J/4>, K xi K 4>'> vj/iPN = aj/Tj,N,a x N,Vil>'N and M = Mj^, M x ,M r 
respectively. The remaining problem is whether formula @ is valid for \ and -0'. There are data for {p'j) I ^i' ph 
and they shows that {p^p)^ pb ~ lA(p? r ) p, / Pb ■ So, we assume that the above is also true for (p^)^/ S (e), i.e. 

<J&)$ 3 (e) = lA<&)jf v (e) (24) 
with (Pr)j^(e) given by (||) . For \ we believe that the inequality 

<&)1fc < (Vt)x B < (Pt)$? (25) 

should be valid and therefore assume that ( |24| ) is true also in this case. Anyway, the exact form of (p^)^ s (e) 
or (Px)^, P ( e ) is n °t ver y important because we checked that the suppression depends on this form very weakly. 
First, we put Kj/^ — and the resulting J/^ff survival factor (for direct J/vP's) differs only a few percent for 
the highest eo from the one calculated with formula (H) unchanged. Second, when we use expression (p[) also 
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for x and ip' ', the evaluated suppression factor is the same as that calculated with the use of (^4|), as far as plots 
are concerned. 

To complete our estimations we need also values of cross-sections for x ~ baryon and tp' — baryon scatterings 
(we will still hold that — meson cross-section is | of x(VO — baryon cross-section). Since J/\I/ is smaller 

than x or tp', x~ baryon and ip' — baryon cross-sections should be greater than J "/ W — baryon one. For simplicity, 
we assume that all these cross-sections are equal. This means that we underestimate J/^ suppression, here. 
The final results of calculations of ( |l3| ) are presented in Figs.|||| for various sets of parameters of our model 
(which are Tf. .,n%,<7b). We performed these calculations for two values of Tf. . — 100, 140 MeV which agree 
fairly well with values deduced from hadron yields ||l0|| . For comparison, also the experimental data are shown 
in Figs, pl-pl The experimental survival factor is defined as 



AC 



—txb — 







(26) 



where 



w _jj± 

AB{ PV ) 



is the ratio of the J/*5> to the Drell-Yan production cross-section in A-B(p-p) interactions 

times the branching ratio of the J/vE' into a muon pair. The values of the ratio for p-p, S-U and Pb-Pb are 
taken from li^B . Note that since the equality o^y = ° V dy ' has been confirmed experimentally up to now 
formula (|26f) reduces to 



AB 

AT - aj /^ 

JV exp - ABaPP 



(27) 



DY 



which is also given as the experimental survival factor, for instance, in |40 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 
£ [GeV/fm 3 ] 
FIG. 5. J/9 suppression in the longitudinally ex- 
panding hadron gas with the " infinite" transverse size 
and for n% = 0.25 fin" 3 and T f . a . = 140 MeV: (a) 
a b — 3 mb, a m = 2 mb; (b) a b = 4 mb, cr m = 2.66 
mb; (c) a b = 5 mb, a m = 3.33 mb; (d) a b = 6 mb, 
<r m = 4 mb . The black squares correspond to the 
NA38 S-U data [^3|, the black triangles correspond 
to the 1996 NA50 Pb-Pb data, the white squares to 
the 1996 analysis with minimum bias and the black 
points to the 1998 analysis with minimum bias ]l3 |. 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 
£ [GeV/fm 3 ] 
FIG. 6. Same as Fig.|E| (except case (d), which is 
not presented here) but for Tf. . = 100 MeV. 



Coming back to examination of our first results presented in Figs.gJ- 
patterns of suppression which agree with the experimental data fairly we 



q, we can see that we have obtained 
1 for some values of parameters of our 



model. The data prefer = 5 — 6 mb and (or) Tf closer to 100 MeV. Note that the dependence on the 
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initial baryon number density is substantial but for higher values of n° B , rather. The lower the initial baryon 
number density, the deeper the suppression. There are two reasons for such a behaviour: the first, for the 
higher baryon number density, there are less non-strange heavier mesons p, u> in the hadron gas of the same 
eo, but these particles create the most weighty fraction of scatters, for which reaction |l]) have no threshold at 
all; the second, the freeze-out time i/. . decreases with increasing n° B for a given eq in our model. For instance, 
for e = 3.5 GeV/fm 3 and T f . . = 140 MeV we have a = 0.172, 0.175, 0.183 and t f . . = 15.7, 14.7, 11.6 fm 
for n° B = 0.05, 0.25, 0.65 fm -3 , respectively. We can see also that the value Ob = 3 mb is too small to obtain 
results comparable with the data, so we will leave aside this value in further investigations. 




0.0 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

£ [GeV/fm 3 ] 
FIG. 7. J/^l suppression in the longitudinally ex- 
panding hadron gas with the "infinite" transverse 
size and for cr& = 4 mb, <j m = 2.66 mb and 
T } . . = 140 MeV. The curves correspond to n° B = 0.05 
fm -3 (solid), n% = 0.25 fm" 3 (short-dashed) and 
n° B = 0.65 fm -3 (dashed) The black squares repre- 
sent the NA38 S-U data [ J13| , the black triangles rep- 
resent the 1996 NA50 Pb-Pb data, the white squares 
the 1996 analysis with minimum bias and the black 
points the 1998 analysis with minimum bias JL4). 




e [GeV/fm 3 ] 



FIG. 8. Same as Fig.| but for T f . . = 100 MeV. 



Now we will include the finite-size effects into our model, i.e. we will take into account that the realistic 
hadron gas has a finite transverse size. This will be done in form of the rarefaction wave moving inward S e t / 
with the sound velocity c s . How to obtain this velocity has been mentioned in Sec. 2 (see also p2|) . With the 
finite-size effects included, the final expression for J/ty survival factor Afh.g.( e o) wu l be given by (|15| ) . To make 
our investigations much more realistic we will also include the possible J/\f' disintegration in nuclear matter, 
which should increase J/\& suppression by about 10% fll5| . But to draw also S-U data in figures, instead of 
multiplying ftfh.g. by N n . m . given by (|7|), we divide N exp by appropriate Mn. m ., he. we define "the experimental 
J/\l/ hadron gas survival factor" as 

Nexp = exp {aj/^ N p L} ■ Nexp ■ (28) 

as the experimental data. 



and values of this factor are drawn in Figs.|9|-|10] and Figs.|12j-|13 

We shall consider the Woods-Saxon nuclear matter density distribution ]36|] , here. The results of numerical 
estimations of (|l5| ) and (28) are depicted in Figs. |§-[l0] for two values of the charmonium-baryon cross-section 



(Jb = 4, 5 mb and the initial baryon number density n B = 0.25, 0.65 fm -3 . The curves for n B = 0.05 fm -3 
almost cover the curves for n B = 0.25 fm~ 3 , so for clearness of the figure we do not draw them. The two values 
of the speed of sound are the maximal values of this quantity possible in the range [7/. . = 140 MeV, To jmax ] 
for the above-mentioned two cases of n B . In fact, we have checked that the results almost do not depend on c s 
allowed in the range. 

It has turned out also that in the case of the transverse expansion, the results almost do not depend on the 
Tf.o. (f° r Tf. . € [100, 140] MeV). This is because the freeze-out time resulting from the transverse expansion, 
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tf.o., trans = Ra/cs (if we assume a central collision and c s constant), is of the order of the freeze-out time 
resulting from the longitudinal expansion for X/. . = 140 MeV. Namely, for Pb and c s = 0.45 we ha,vetf. . t t ra ns — 
15.8 fm, which is very similar to values of £/. . for 7y o . = 140 MeV given earlier. For 7y o . = 100 MeV, 
i/.o. = 111.0, 101.0, 72.5 fm for n% = 0.05, 0.25, 0.65 fm -3 respectively, so the hadron gas ceases because of 
the transverse expansion much earlier. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 
£ [GeV/fm 3 ] 
FIG. 9. J/\& suppression in the longitudinally 
and transversely expanding hadron gas for the 
Woods-Saxon nuclear matter density distribution and 
a b = 4 mb, a m = 2.66 mb and T f . . = 140 MeV. The 
curves correspond to n° B = 0.25 fm~ 3 , c B — 0.45, 
r = 1.2 fm (solid), n% = 0.65 fm" 3 , c s = 0.46, 
r = 1.2 fm (dashed) and n% = 0.25 fm -3 , c s = 0.45, 
ro = 1.12 fm (short-dashedh The black squares rep- 
resent the NA38 S-U data fllJ , the black triangles rep- 
resent the 1996 NA50 Pb-Pb data, the white squares 
the 1996 analysis with minimum bias and the black 
points the 1998 analysis with minimum bias [Q , but 
the data are "cleaned out" from the contribution of 
J/ty scattering in the nuclear matter in accordance 
with (||). 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 
e [GeV/fm 3 ] 
FIG. 10. Same as Fig. ^ but for at = 5 mb and 
m — 3.33 mb. 



Generally, taking into account also the transverse expansion changes the final (theoretical) pattern of J/^ 
suppression qualitatively. First of all, the curves for the case including the transverse expansion are not convex, 
in opposite to the case with the longitudinal expansion only, where the curves are. But still, theoretical curves 
do not fall steep enough at high e to cover the data area completely. Nevertheless, from Figs. p|-|To| we can see 
that for some choice of parameters, namely for <t& somewhere between 4 and 5 mb and for ro = 1.12 fm, we 
would obtain a quite satisfactory curve. And we should remember that since we have one overall charmonium- 
baryon cross-section ab, our final results underestimate the suppression (for \~ > "0' — baryon scattering the 
cross-section should be greater than for J/W). To support our conclusion in more visible way we present main 

results from Figs.|^-[l0] in Fig.jll], where original data H] for and J/^ survival factors given by (|^) 

jr-^ and now as functions of Et are presented. The change of the variable from eg to Et has 



multiplied by 

been done with the use of (|l9|) and 6q = eo(b) expressed by the solid line in Fig. 
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50 t 
45: 




E T [GeV] 



FIG. 11. J/SP survival factor times pp'^ in the longitudinally and transversely expanding hadron gas for the 

<T r>i' 

Woods-Saxon nuclear matter density distribution and n B = 0.25 fm~ 3 , Tf. . = 140 MeV, c s — 0.45 and ro = 1.2 fm. The 
curves correspond to ab — 4 mb (solid) and Ob = 5 mb (dashed). The black triangles represent the 1996 NA50 Pb-Pb 
data, the white squares the 1996 analysis with minimum bias and the black points the 1998 analysis with minimum bias 




1.2 - 




0.0 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

e [GeV/fm 3 ] 
FIG. 12. J /SP suppression in the longitudinally 
and transversely expanding hadron gas for the 
Woods-Saxon nuclear matter density distribution and 
a b = 4mb, a m = 2.66 mb, T f . . = 140 MeV and 
ro = 1.2 fm but for eo(6) given by ([ill). The curves 
correspond to n° B = 0.25 fm -3 , c s = 0.45 (solid) and 
n° B = 0.65 fm" 3 , c s = 0.46 (dashed^ The black 
squares represent the NA38 S-U data 1l3j , the black 
triangles represent the 1996 NA50 Pb-Pb data, the 
white squares the 1996 analysis with minimum bias 
and the black points the 1998 analysis with minimum 
bias [Q, but the data are "cleaned out" from the con- 
tribution of J/S& scattering in the nuclear matter in 
accordance with (Ea). 
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FIG. 13. Same as Fig.|l^ but for crj, = 5 mb and 

o m = 3.33 mb. 



Note that the main disagreement with the data reveals in the last experimental point of the 1998 analysis 
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Q (see Figs.|-0). But the error bar of this point is very wide. Additionally, there is some contradiction in 
positions of the last three points of the 1998 data. In fact, the middle point goes up above the left one and the 
last falls far below the two earlier. This suggests that more experimental data is needed in the region of high 
eo to state definitely whether the abrupt fall of the experimental suppression factor takes place or not there. 

Now we repeat again the numerical estimations of (|l^) for the case with the finite-size effects and the Woods- 
Saxon nuclear matter density distribution included, but for the dependence of eo on b given by (pl|). The results 
are presented in Figs.[l2T3. Note that the theoretical curves are two-valued around eo = 2.8 GeV/fm 3 . This 
is the result of our approximation of eo(b) by (pT|). This expression allows for two different values of b, which 
give the same eo in some range of the impact parameter b. This is shown in Fig. |] (short-dashed line) . We can 
see that for b < 7.9 fm there are two different values b\ and &2 such that eo(£>i) = epfe). This causes that 
results plotted in Figs. p~2| - [T^ are qualitatively different from those presented in Figs.||-[lO] in the region close to 
the maximal eo reached at the collision. Comparing Figs.p|-^0| with Figs.[l2]-|l3|, we can see also that the pattern 
of J/vP suppression depends on the shape of eo as a function of b. Clarification of this dependence would be 
very helpful to obtain more realistic picture of the J /^Sf dissociation in hadron medium during the heavy-ion 
collisions. 

As a final remark, we think that it is difficult to exclude J/vf scattering in the hot hadron gas entirely, as the 
reason for the observed J/^S suppression at this point (see also In our model the most crucial parameter is 

the charmonium—baryon inelastic cross-section and the final results depend on its value substantially. Therefore 
it is of the greatest importance to establish how this cross-section behaves in the hot hadron environment. Some 
work has been done into this direction p2-H3|], but results presented there differ from each other and are based 
on different models. However, the newest estimations of tt + J/^, p + J/^ and J/*f> + N cross-sections at high 
invariant collision energies |^5|,^6| agree with the values of o~b and a m assumed in our model. We would like to 
add also at this point that the charmonium — hadron inelastic cross-sections have been considered as constant 
quantities here. For sure, they should not be constant and the results of just mentioned papers suggests 
that they are not, indeed. They are growing functions of the invariant collision energy y/s. So, the naive 
reasoning should direct us to the conclusion that the increase of eo (or in other words Et) causes the increase of 
the invariant collision energy yfs on the average and further the increase of the charmonium — hadron inelastic 
cross-sections. This could influence the final patterns of J/^f suppression in such a way that J survival factor 
would behave according to the solid curve of Fig.pl] for low eo (Et) but then, as the charmonium — hadron 
inelastic cross-sections would increase with eo {Et), the factor would go closer to the dashed curve of Fig.0 
for high e (Et)- So, the experimental pattern of J suppression could be recovered in this way. Therefore, as 
a final conclusion we can say that it is difficult to ruled out the conventional explanations of J/^f suppression 
completely, at present. 

We would like to stress again that the behaviour of the experimental J/\P suppression factor at high Et (or 
otherwise at high eo) has not been clear yet. In fact, the abrupt fall of this factor (what could suggest the 
appearance of the quark-gluon plasma) is indicated only by the one point (the last) of the 1998 NA50 analysis 
fl4|j . So, to draw a definite conclusion more experimental data far above Et — 120 GeV are needed. This region 
will be reached in upcoming RHIC runs and their results should answer the question: is the J/\& suppression a 
signature of the existence of the quark-gluon plasma, or not?. 

Note added. When our paper was completed we became aware of [E7| where the twin figure (denoted as Fig. 5 
there) to our Fig.^ was presented. But the appearance of the quark-gluon plasma is the main reason for J/ty 
suppression there. It is also claimed that results shown in that figure "provide evidence for the production of 
the quark-gluon plasma in central high-energy Pb-Pb collisions" . This entirely confirms our conclusion that the 
status of J/^S suppression as a signal for the quark-gluon plasma appearance is far from being clear at present. 
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